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Abstract — We consider a communication system where two 
transmitters wisli to excliange information ttirougti a Iialf-duplex 
relay in the middle. The channels between the transmitters and 
the relay have asymmetric channel gains. More specifically, the 
channels are assumed to be synchronized with complex inputs 
and complex fading coefficients with an average power constraint 
on the inputs to the channels. The noise at the receivers have 
the same power spectral density and are assumed to be white 
and Gaussian. We restrict our attention to transmission schemes 
where information from the two nodes are simultaneously sent to 
the relay during a medium access phase followed by a broadcast 
phase where the relay broadcasts information to both the nodes. 
An upper bound on the capacity for the two phase protocol 
under a sum power constraint on the transmit power from all 
the nodes is obtained as a solution to a convex optimization 
problem. We show that a scheme using channel inversion with 
lattice decoding can obtain a rate a small constant 0.09 bits from 
the upper bound at high signal-to-noise ratios. Numerical results 
show that the proposed scheme can perform very close to the 
upper bound. 

I. Introduction 

In our previous work [3], we studied the bi-directional relay 
problem where two users exchange information with each 
other through a relay in the middle. Therein we assumed 
same channel gains between the nodes and the relay. In this 
work we look at a more practical scenario with asymmetric 
channel gains. There has been recent work on the bi-directional 
relay problem with asymmetric channel gains. In [1], sev- 
eral schemes including compress and forward, amplify and 
forward, decode and forward and mixed forward have been 
suggested. In our work, we show the benefit of using nested 
lattices for this problem. There has been some work on using 
nested lattices for asymmetric channel gains in [5,6]. In [6], 
a scheme that is optimal at high signal-to-noise ratios (snr) is 
given for a given realization of the asymmetric channel. 

In this paper, we consider a scenario of a fading channel 
with L channel realizations and assume that channel gains are 
known to all the nodes. The transmit power is adapted as a 
function of the channel gains. We show that an upper bound 
to the achievable exchange capacity (defined later) can be 
obtained as the solution to a convex optimization problem. We 
also show that symmetric nested lattices with an appropriate 
power allocation policy can perform close to the upper bound 
at high SNRs. [] 

This work was supported by the National Science Foundation under grant 
CCR-05 15296. 



II. System Model and Problem Statement 




Fig. 1. Problem setup with fading links ha and hf, 

We study a simple 3-node linear Gaussian network, with 
asymmetric channel gains as shown in Fig. [T] Node A and 
node B wish to exchange information between each other 
through the relay node R. However, nodes A and B cannot 
communicate with each other directly. The nodes are assumed 
to be half-duplex, i.e. a node can not transmit and listen at 
the same time. We use the block fading model to model 
the channels between the nodes and the relay. The channel 
gains between the nodes and the relays remain constant over 
a coherence time interval. It is assumed that the transmis- 
sion happens over L such coherence time intervals. Each 
coherence interval corresponds to N uses of the channel. 
Hence, effectively NL uses of the channel are available for 
communication. For the channel between the node A and 
the relay R, each of the coefficients of the L length vector 
G C^, represent the channel gain for each coherence time 
interval. Similarly the channel between the node B and the 
relay R, the relay R and node A and the relay R and node 
B are captured by the channel gain vectors hbr, hra and 
hrb G respectively (vectors are denoted by bold face letters 
such as V throughout the paper). The channels are assumed to 
be reciprocal, i.e., har ~ hra ha and hbr = hrb := hb. 

Let e {0, l}'^'"^ and ub <E {0, l}'-''-^ be the information 
vectors at nodes A and B. We assume a protocol where the 
communication takes place in two phases at each coherence 
time interval i { i E {1,2 . . . L}). The phases are the multiple 
access (MAC) phase and the broadcast phase. A G [0, 1] is the 
fraction of channel uses for which the MAC phase is used and 
(1 — A) is the fraction of channel uses for which the broadcast 
phase is used. It is assumed that communication in the MAC 



• A scheme is proposed which at high snrs is away from 
the upper bound by at most 0.09 bits(see Theorem 3). 
The scheme uses nested lattice encoding and the transmit 
power is chosen to be inversely proportional to the 
channel gains. 

IV. Upper Bound for the two phase protocol 

We can easily obtain an upper bound for the two phase 
protocol using cut-set arguments and as in [1]. In our problem 
model the channel remains constant in each coherence time 
interval. Hence, the channel over L such coherence time 
intervals can be modeled as a set of L parallel channels. At 
each coherence time interval i, the maximum information rate 
that can be transmitted from node A to node B is bounded by 
the minimum of the information capacity from node A to relay 
node R, and relay node R to node B. This can be expressed as, 
min{AC(|/iai|2p„i),(l - A)C{\hbi\''Pri)}, where C{x) := 
log(l + x). Here A is the fraction of time node A transmits 
and (1 — A) is the fraction of time the relay node transmits. 
Similarly the rate that can be transmitted from node B to node 
A is bounded by the minimum of the information capacity 
from node B to relay R and from relay R to node A, which can 
be expressed as min{AC(|/i(,ipP6j), (1 — A)C(|/ioiP-Pri)}- 
Hence the total rate that can be transmitted over L such 
coherence time intervals can be expressed as the sum of the 
rates at each coherence time interval. Combining the above, 
the upper bound on the exchange capacity can be expressed 
as the solution of an optimization problem given by. 



and broadcast phases are orthogonal. For example, this could 
be in two separate frequency bands (or in two different time 
slots) and hence the MAC phase and broadcast phase do not 
interfere with each other. 

a) MAC phase: During the MAC phase of each coher- 
ence time interval i, nodes A and B transmit while the relay 
node listens. Xai G and Xbi € are the transmitted 
vectors at nodes A and B, respectively. The MAC phase takes 
place in I^N uses of the complex additive white Gaussian 
Channel (AWGN) channel. Further it is assumed that the two 
transmissions are perfectly synchronized. Hence the received 
signal at the relay yri G C^^, is given by 

where the components of Zri G are independent 

identically distributed (i.i.d) complex, circularly symmetric 
Gaussian random variables with zero mean and unit variance. 
The average transmit power at the nodes A and B in the z*** 
coherence time interval is given as = Pai and 

-^[ll^bill ] = Pbi- 

b) Broadcast phase: During the broadcast phase, the 
relay node transmits x^i G C'^~^^^ in the i*'' coherence time 
interval to both nodes A and B. The nodes A and B receive 
yai and ybi, respectively where 

Ybi = ^fejXbi -I- Zbi 

The average transmit power at the relay node during the i*'' 
coherence time interval is given by Pru and the receiver noise 
at the two nodes is complex Gaussian with zero mean and unit 
variance. 

Further, it is assumed that there is a total sum power 
constraint over all the nodes. Since the MAC phase is used 
during the fraction A and the broadcast phase is used during 
the fraction (1 — A) of the available time slots, the total power 
constraint is expressed as 

Ax^ A^ (1-A)^ 

i—l i — 1 i—1 

We are interested in power allocation strategies and good 
encoding/decoding schemes that maximize the amount of 
information (maximize k,J{LN) + ki,/{LN)) that can be 
exchanged rehably (such that the probabihty of error can be 
made arbitrarily small in the limit of TV — » c»). We refer to the 
maximum value of ka/{LN) + ki,/{LN) that can be reliably 
exchanged with a given scheme as the exchange rate for that 
scheme. The exchange capacity is then the supremum of all 
such rates over the encoding schemes. 

III. Main Results and comments 

We only consider the case when A = 0.5, or the MAC and 
broadcast phase each use the channel half the time. For this 
case, the main results in this paper are 

• An upper bound on the exchange rate is setup as a convex 
optimization problem. 



1 

maximize - V min{AC(|/iai|2p„i), (1 - A)C{\hM\''Pri)} 
^ i=i 

1 ^ 

+ j^^mm{AC{\hMfPM), (1 - A)C{\KifPri)} 

i=l 

subject to ^E-^«+xE'^*'' + ^^;^E^"-'^ 

Pai, Pbi,Pri>0,ie 1,2... L. (1) 



Since A is fixed, we can see that the objective function given 
above is concave and also the constraints form a convex set. 
Hence convex optimization techniques can be easily applied 
to this setup to get the optimal solution for the above convex 
problem. 

Though the problem is convex, it is tough to get analytical 
results for general snr. However to develop some intuition 
and for ease of analysis, we obtain the power allocation 
strategy by solving the optimization problem under the high 
snr approximation, i.e. C{x) := log(a;). We make this precise 
in the following theorem. 

Theorem 1: For finite L and A = 0.5, and under the high 
snr approximation, the optimal power allocation for the upper 
bound on the achievable rate for different ranges of /t| := 
|/iaiP/|/i6iP is given by. 




Fig. 2. Power allocation as a function of k for the upper bound 



Case 1: < Kf <= ^i±^ , 
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Case 3: > 



Pa, = Pz 
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Phi ~ P Pri P - 



Proof: The problem can be solved analytically using the 
method of Lagrange multipliers and applying the Karush-Kuhn 
Tucker conditions as given below. The optimization problem 
for finite L and A = 0.5 can be expressed below as, 



maximize 



subject to 



\c{\hu?Pu) > Rtl \c{\Kr?Pri) > R^ 



1 ^ 



Pa^,PM,Pr^ > 0, Z G { 1 , 2 . . . L} 



We next use the high snr approximation to simplify the 
analysis. Under the high snr approximation C{x) w logx. 



The Lagrangian can then be expressed as, 



C = -- 



2 = 1 

{Pa^ + Pm + Pr^)- 2P^ 

L 

Y,p^{2B!{l-\og{\K,^P^ 

L 

Y.y^f{2R^2-\og{\hu\''Pr.)) 



L 



5]4*)(2i?«-iog(|/i,,pn,)) 



i=l 
L 



+ 5^MW(242-log(|/i„pP„)) 



i=l 



Next taking the partial derivative of C with respect to each 
variable in £ and equating to zero, gives us the following set 
of equations. 



1 
L 



2fi['^ + 2nf = 0, - - + 2/4'^ + 2^, 



1 



,( ■0 - 



L 







A 



^h_ ^ _ = — - 

L Pai L Pfji L 

Pr 



= 



Solving for A from the above set of equations gives A = ^p- 
Solving again the above set of equations for Pai, Pm and P^i 
and together with the Karush-Kuhn Tucker(KKT) conditions, 
gives us the required power allocation as a function of k^. ■ 

In the above proof, the parameter A is not a function of 
the channel gains. This is so, since we have made the high 
snr approximation. This implies that the total power allocated 
Pai + Pm + Pri, during each coherence time interval remains 
the same. However, the power in the individual nodes will 
vary based on the channel gains. From Theorem 1, we can 
see that in case 1 we have at low nf, n^Pai ~ Pm- In case 3, 
where k? >> 1, nfPai ~ Pm- Also for case 2, where k| w 1, 
lifPai ~ Pm- In other words this implies that \hai\'^Pai ~ 
{hul^PM- In the next section, we propose a scheme that makes 
use of this property to obtain results close to the upper bound. 

V. Achievable Scheme using channel inversion 

AND LATTICE CODING 

In this section we discuss our achievable scheme based on 
nested lattice decoding by Erez and Zamir[2]. The proposed 
scheme follows closely the lattice coding scheme discussed in 
our previous work [3] and in [4]. The main idea in [3], is that 
suppose each of the nodes A, B and the relay R, has the same 
power constraint (say Pa), then at high signal to noise ratios, 
a rate close to the upper bound C(Pa) can be exchanged. 
In other words, the nested lattice coding approach works best 
when the receiver channel signal strengths from the two nodes 
are the same. 



In our proposed scheme, we enforce \hai\'^Pai — l^biP-Pfci 
for every i*^ coherence time interval in the MAC phase. This 



matches closely with the observation in section IV where 
the power allocation profile at the nodes satisfies \hai\'^Pai ~ 
\hbi\'^Pbi- This means that each node uses a coarse lattice of 
power P\i, and each node performs a channel inversion at the 
transmitter, so that the relay receives equal signal strengths 
from both the nodes. For the broadcast phase, to ensure that 
the nodes A and B can decode the message from the relay, 
we enforce that transmit power at the relay is always larger 
than the transmit power at the nodes, or P^j > Paii Pu- 

First let us explain our achievable scheme in detail for the 
i*'' coherence interval with known channel gains hai, fiM- 
We first obtain the power allocation profiles Pai, Pm and Pri, 
based on the additional requirements of \hai\'^Pai — \hbi\'^Pbi 
and Pri > Pai, Pbi- The allocation is discussed in more detail 
later in the section. Let us next define P^i := \hai\'^Pai — 
\hbi\'^Pbi and for each coherence time interval i, choose a 
nested lattice structure having a fine lattice A{ with a coarse 
lattice Ai nested in it. The second moment per unit dimen- 
sion of the coarse lattice is PAi/2. Also the channel model 
considered in this problem setup has complex inputs and 
complex noise, when compared to the real Gaussian channel 
model in [3]. The complex channel provides two degrees of 
freedom. To take advantage of this we can perform nested 
lattice coding separately along the in-phase and the quadrature 
phase components. In all the vectors discussed below namely 
Xai, Xbi, tai, tbi, Uai and Ubi are complex vectors and can be 
expressed as the complex sum of their in-phase and quadrature 
phase components. For example Xai can be expressed as 
Re{xai} + jIm{Xai}. 

First at each coherence interval i during the MAC phase, the 
data at the nodes A and B are mapped to lattice points tai ™d 
tbi respectively. Let Uai and Ubi be dithers that are uniformly 
distributed over the coarse lattice. The in-phase component of 
the dither and the quadrature phase component are independent 
of each other and each distributed uniformly in the coarse 
lattice Ai, with second moment P\i. An output (tai — Uai) 
mod Ai is obtained at node A and (tbi— Ubi) mod A^ at node 
B. Here ti mod A^ represents (Re{ti} mod Ai)+j(Im{ti} 
mod Ai) 

Hence the transmitted vector Xai at node A is given by, 
(tai - Uai) mod Ai 



(2) 



Here the numerator (tai — Uai) mod Ai is the sum of 
the in-phase and quadrature phase components, expressed as 
Re(tai — Uai) mod Ai + jlm(tai — Uai) mod Ai. Hence 
the second moment of the numerator per unit dimension is 
PAi/2 + PAi/2 = Paj- Hence the average transmit power on 
Xai is PAi/|/iaiP — Pai- Thus wc meet the average power 
constraint of Pai at node A. Similarly the transmitted vector 
Xbi at node B is given below. This also meets the average 
power constraint Pbi- 

(tbi - Ubi) mod Ai 



The relay receives 

Yri = /lajXai + /ifaiXbi + Zri (4) 

or 

Yri = (tai-Uai) mod Aj + (tbi - Ubi) modAi+Zri(5) 

The decoder next forms (yri + Uai + Ubi) mod Ai and per- 
forms nested lattice decoding and decodes to tri = (tai + tbi) 
mod Ai with high probability, as long as the transmission rate 
from each of the nodes is less than 2{i log(0.5 + PAi)}- The 
factor 2 is present because the channel coefficients are complex 
and we have 2 degrees of freedom. 

The relay next forms (tri — Uri) mod Ai and during the 
broadcast phase transmits 



P'ri 



{(tri -Uri) mod Ai} 



The relays can decode to tri as \hai\'^Pri,\hbi\'^Pr 



(6) 

> Pa,, 



since Pri > Pai, Pbi- Hence effectively a rate of ilog(0.5 
PAi) can be achieved by the nodes. 

Also define D{x) :— w.c.e{log(0.5 + a;),0.51og(l + 
2x)},x > 0. Here u.c.e denotes the upper concave envelope 
of the two functions. Hence the optimization problem can be 
expressed as follows with a few more constraints added. 



maximize 



1 

- mm{AD{\har\^Pa^),{l - A)D{\hb^\^ Pr^)} 



^Y.mm{AD{\hb4'Pb^), (1 - A)i?(|/i,,|^P„)} 

1=1 



subject to 



Pa 



It 

1=1 
\h PP 

P > P 

Pri ^ -Pbi 5 

Pa^,Pb^,Pr^>0,iel,2. 



(1-A) 



i=l 
— \hbi\'^Pbi, 



J2P"<P, (7) 

i=l 

(8) 
(9) 

(10) 



.L 



The above optimization problem is solved for the case 




Xbi 



(3) 



Fig. 3. Power allocation as a function of k for our proposed scheme 

for A = 0.5 under the high snr assumption with D{x) 
approximated by log(a;). The next theorem gives us the power 



allocation profile at the different nodes for different values of 

:= \K^\V\hb^?■ 

Theorem 2: For finite L with A ~ 0.5 and under the 
high snr approximation, the optimal power allocation for the 
achievable scheme for different ranges of k| := |ft.aiP/|^fciP 
is given by. 

Case 1: Q < nj <= 1, 

2 2 

p . — p r p. . — p ' p . — p r 



Cflie 2; > 1, 



P- 



1 + 2k,^ 



2At? 



2k? 



1 + 2k,^ 



1 + 2k? 

Proof: The problem can be solved analytically using the 
method of Lagrange multipliers and using the Karush-Kuhn 
Tucker conditions following along the same lines as in the 
proof of Theorem 1 . ■ 

VI. Comparison of the upper bound and the 

ACHIEVABLE SCHEMES 

Next we state the main theorem of this paper that compares 
the upper bound and the achievable scheme for known channel 
state information in L coherence intervals. 




- Upper Bpund 

- Channel inversion sciieme with Lattice cotling 

- Amplify aniJ ForwartJ scheme 



Transmit Power P 



Fig. 4. Comparison of bounds for A = 0.5 

Theorem 3: For general L with A = 0.5, and under the 
high snr approximation, the achievable rate using the channel 
inversion scheme with lattice decoding suffers at most a 
constant 77 = 0.08972 bits per complex channel use from the 
upper bound. 

Proof: We compare analytically the results of Theorem 1 
and Theorem 2 under the high snr approximation. For different 
values of Ki we compare the exchange rates and we can easily 
show that the channel inversion scheme suffers at max 0.08972 
bits per complex channel use from the upper bound. ■ 
Theorem 3, hence captures the loss due to adding the 
additional constraints to the upper bound, and the loss is 
found to be really small. Fig |4] shows some results obtained by 
numerical solving the optimization problems without the high 
snr approximation. Here we have L = 100, and the channel 
coefficients are taken from a Rayleigh distribution, and they 
have unit variance. The channel coefficients are fixed and we 
evaluate the average rate per channel use, each for the upper 
bound, the lattice based scheme and also the amplify forward 
scheme. 



Vll. Practical issues with power allocation 

DESIGN 

In the previous sections, we discussed techniques to com- 
pute the optimal power allocation for a given ha, hb and 
hr. However, in practice we are interested in maximizing the 
average exchange capacity, i.e., the problem is to 

maximize Eh^^, [mm{AC {\haf Pa) , (1 - A)Ci\hbfPr)} 

+ mm{AC{\h\^Pt),{l-A)C{\ha\^Pr)}] 
subject to Eh^j,, [APa + APb + (1 - A)P,] < P 

Pa,Pb,Pr > 0. 

Due to the ergodic nature of the channel, the above optimiza- 
tion problem is identical to the one in ([T]) when L ^ 00. The 
result in Fig.|4]have been obtained by solving the optimization 
problem for one realization of ha,hb and hr but with L = 
100. Note that while computing the optimal power allocation 
policy requires us to use a large value of L and optimize the 
policy, once this policy is fixed, the actual transmit power is 
chosen based only on the instantaneous channel realization. 

Vlll. Conclusion 

In this paper we studied the bi-directional relay problem 
in which the channels between the nodes and the relays 
were assumed to have complex inputs with complex fading 
coefficients. We studied power allocation policies at the nodes 
for a two phase transmission scheme under the sum transmit 
power constraint over all nodes. For A = 0.5, where each 
phase uses the channel exactly half the time, we obtained 
an upper bound on the exchange capacity as a solution to a 
convex optimization problem. We proposed a scheme using 
nested lattice encoding with the transmit power chosen to 
be inversely proportional to the channel gains. We obtained 
analytical solutions for the exchange capacity under the high 
snr approximation and showed that our proposed scheme can 
obtain a rate which is at most 0.09 bits away from the upper 
bound. For A 7^ 0.5, we were unable to obtain a good 
performance using a simple channel inversion power allocation 
policy. However, it can be shown that using lattice codes with 
asymmetric rates [6] at the nodes, the upper bound can be 
achieved at high snrs. 
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